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The  long  range  objectives  of  this  research  are  to  identify  and  correlate 
important  properties  of  the  heat  transfer  process  that  can  be  used  to  optimize  the 
heat  sink  capabilities  of  evaporating  thin  films.  Initially,  we  proposed  to  emphasize 
the  "Interline  Heat  Flow  Parameter".  Using  the  results  of  our  theoretical  research  we 
now  plan  to  enhance  this  objective  and  include  the  evaluation  of  the  "Capillary 
Pumping  Parameter"  for  completeness.  The  interline  heat  flow  parameter  is 
important  in  ultra-thin  films  where  the  liquid-solid  interfacial  force  field 
predominates,  whereas  the  capillary  pumping  parameter  is  important  in  thicker  films. 

(_ 

yAh  \  ft  '  AAh 
v  )  3  ^  v 

where 

2rrA  is  the  Hamaker  constant, 

ft  is  the  thickness  of  the  liquid  film. 

Ah  is  the  heat  of  vaporization, 

y  is  the  surface  tension, 

K'  is  the  curvature  gradient, 

v  is  the  kinematic  viscosity, 
yAh 

'  is  the  capillary  pumping  parameter, 

AAh 

'  is  the  interline  heat  flow  parameter, 
with  K1  and  ft’  evaluated  at  8(x). 

Therefore,  the  experimentally  measured  film  thickness  profile  can  be  used  to 
calculate  the  integral  heat  sink. 


There  are  two  main  aspects  of  the  program:  experimental  and  theoretical.  We 
will  discuss  the  experimental  portion  first  by  outlining  significant  decisions  and 
accomplishments  that  have  occurred  during  the  first  16  months.  This  will  be  followed 
by  a  discussion  of  the  theoretical  accomplishments. 

EXPERIMENTAL:  GENERAL 

The  major  experimental  task  completed  during  this  period  was  the  design  and 
fabrication  of  a  heat  transfer  cell  for  studying  thin  film  heat  transfer  and  fluid  flow  in 
the  contact  line  region.  A  considerable  amount  of  effort  was  directed  at  designing  and 
testing  a  new  capillary  feeder  for  the  heat  transfer  cell.  This  concept  was  tested  with 
various  flow  geometries  and  it  was  found  that  a  circular  capillary  feeder  works  best 
in  terms  of  flow  symmetry,  absence  of  "edge  effects",  efficient  sealing,  and 
symmetrical  temperature  and  heat  flux  profiles.  We  have  decided  to  evaluate  the 
integral  evaporation  rate  in  the  contact  line  region  by  measuring  the  total  fluid  flow 
rate  from  the  capillary  feeder  and  evaluating  this  in  light  of  the  measured  thickness 
profile  using  Eq.  (1).  bins  eliminates  recondensation  uncertainty.  Fig.  1  is  a  schematic 
view  of  the  circular  capillary  feeder  without  its  enclosure.  A  photograph  of  the 
complete  cell  is  given  in  Fig.  2. 

Both  the  ultra  thin  (ft  <  20  nm)  and  the  thin  film  (ft  >  20  nm)  regions  of  an 
evaporating  liquid  film  on  a  standard  surface  (bare  Si  wafer)  can  be  studied  using  two 
different  optical  techniques,  elfipsometry  and  interferometry.  The  advantage  of  this 
particular  cell  is  its  flexibility  -  both  the  ellipsometer  and  the  high  power 
interferometer  can  be  used  separately  or,  for  some  magnifications,  simultaneously 
and  specific  design  features  allow  high  heat  flux  studies. 

The  cell  is  basically  a  miniature  capillary  feeder,  the  liquid  spreads  towards  the 
hot  spot,  gets  evaporated  and  is  replenished  from  the  upstream  source  by  the  action 
of  capillary  suction  and  disjoining  pressure.  The  cell  has  the  provision  of  feeding 
liquid  from  outside  and  the  feeding  rate  will  be  monitored  to  evaluate  the  integral 
ipornlion  rate  for  a  particular  power  input.  The  heat  source  is  a  circular  Plati'i'im 
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heater  located  at  the  center  of  the  substrate.  The  temperature  profile  will  be 
measured  by  using  small  thermocouples  attached  at  the  back  of  the  substrate.  We 
have  decided  to  use  thermocouples  instead  of  Silicon  Carbide  thermistor  arrays 
because  a  current  doctoral  thesis  in  our  program  by  E.  Kiewra  demonstrated  that 
there  is  virtually  no  difference  in  the  temperature  data  collected  by  these  two 
different  methods.  Therefore,  the  added  complexity  associated  with  the  thermistors 
is  not  rewarded  by  improved  understanding. 

Another  crucial  factor  in  determining  the  stability  and  equilibrium 
configuration  of  the  liquid  meniscus  on  the  Silicon  surface  is  the  cleanliness  and  liquid 
purity.  The  cleaning  and  the  liquid  purification  procedures  will  be  discussed  in  the 
next  section.  An  important  feature  of  this  cell  is  that  it  can  be  evacuated  to  a  high 
vacuum.  This  will  reduce  the  water  vapor  (which,  if  present  with  organic  test  fluids, 
can  affect  the  interfacial  force  field)  adsorbed  on  the  inner  surface  of  the  cell  walls. 

EXPERIMENTAL:  DESCRIPTION  OF  THE  CELL 

The  substrate  of  the  cell  is  a  single  crystal  Silicon  wafer  7.6  x  10’2  m  in 
diameter  and  about  4  x  1(H  m  thick.  The  Silicon  substrate  has  a  high  degree  of 
surface  uniformity,  high  reflectivity  and  its  optical  and  physical  properties  are  well 
known.  In  the  future  it  will  be  possible  to  use  micromachining  techniques  on  the 
Silicon  wafer  to  make  micro-grooves  to  increase  the  contact  line  length  and  thereby 
enhance  the  heat  transfer  process.  The  heat  source  is  a  circular  shaped  heater  at  the 
center  of  the  wafer  made  by  painting  and  baking  a  Platinum  compound  at  the 
backside  of  the  wafer.  Chromel-Alumel  (type  K)  thermocouples  are  attached  along 
the  centerline  on  the  backside  of  the  substrate. 

Fig.  2  shows  a  photograph  of  the  cell.  The  cell  consists  of  three  parts.  The  top 
part  is  a  hollow  prism  milled  from  a  solid  Aluminum  block.  The  two  sides  are  at  an 
angle  of  70°  with  the  horizontal.  This  is  the  optimum  angle  of  incidence  for 
ellipsometric  measurements.  Each  side  contains  two  viewing  ports  •  supplied  uy  MDC 
Vacuum  Corporation,  made  of  7056  glass  and  they  are  bakeable  to  relieve  residual 
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stresses.  The  viewpoints  are  mounted  on  three  inches  long  nipples  to  keep  them 
away  from  the  evaporation  areas  and  thereby  minimizing  condensation  on  them.  The 
top  contains  another  viewpoint  for  viewing  the  liquid  level  in  the  active  areas  so  that 
the  evaporation  rate  can  be  determined  while  the  process  is  going  on.  It  also  has  a 
connection  for  a  vacuum  line.  Nitrogen  supply  and  a  feeding  system.  The  second  part 
is  an  adjustable  hollow  cylinder  which  fits  the  hole  in  the  top  portion  of  the  prism 
with  matching  threads.  During  interferometric  studies  this  part  can  reach  very  close 
to  the  substrate  surface.  The  microscope  objective  goes  into  the  cylinder.  The 
objective  is  separated  from  the  inside  atmosphere  by  a  heated  glass  plate.  To  allow 
unobstructed  laser  light  to  enter  during  the  ellipsometric  studies  this  assembly  can  be 
lifted  using  the  threads.  This  feature  enables  the  use  of  both  the  ellipsometer  and  the 
high  power  microscope  for  interferometry.  The  third  part  consists  of  an  Aluminum 
plate  which  is  carefully  machined  to  form  the  capillary  groove.  The  groove  holds  and 
feeds  liquid  to  the  circular  meniscus.  The  Silicon  plate  rests  in  between  the 
Aluminum  plate  and  a  Teflon  plate  -  secured  in  place  tightly  by  a  Teflon-coated  O- 
ring  and  a  set  of  12  screws.  Beneath  the  Teflon  plate  is  a  thick  Aluminum  base  plate 
which  holds  the  top  prism  and  the  inside  parts.  A  hole  runs  between  this  plate  and 
the  Teflon  plate  to  reach  the  Silicon  wafer  to  allow  easy  passage  of  thermocouple  and 
heat  source  wires.  The  end  of  the  hole  has  the  provision  for  a  screwed  in  cap  which  is 
kept  closed  during  vacuum  operation  for  cleaning.  The  vacuum  line  and  the  Nitrogen 
line  are  connected  to  the  cell  by  a  three  way  valve.  During  vacuum  operation  the 
Nitrogen  line  is  kept  closed.  Subsequently,  we  close  the  vacuum  line  and  the  Nitrogen 
valve  is  opened  slowly  to  feed  extra-dry  grade  Nitrogen  inside  the  cell  for  operation 
at  a  positive  gauge  pressure.  The  positive  gauge  pressure  minimizes  the  diffusion  of 
impurities  into  the  cell. 

EXPERIMENTAL:  CLEANING 

Cleaning  of  the  cell  assembly  and  particularly  the  Silicon  wafer  is  very  critical  to 
obtaining  reliable  experimental  data.  The  Silicon  substrate  is  very  sensitive  to  its 
cleaning  procedure  and  the  equilibrium  film  thickness  is  a  function  of  its  cleaning 
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history.  For  our  purpose  the  term  "cleanliness"  is  defined  as  the  minimization  of  any 
dust  particles  on  the  Silicon  surface  or  cell  surroundings,  a  well  characterized  Silicon 
surface  (presence  of  oxide,  absence  of  any  organic  layers  etc.)  and  the  absence  of 
anything  which  may  contaminate  the  liquid.  We  will  use  a  standard  cleaning 
procedure  before  each  and  every  experiment  so  that  the  substrate  properties  will 
remain  constant.  We  are  also  evaluating  different  cleaning  processes  and  a  liquid 
purification  system  which  will  be  discussed  later. 

We  will  describe  our  standard  cleaning  process  first.  The  Cell  is  first  rinsed  well 
with  DI  water,  then  immersed  in  pure  Ethanol  and  put  in  an  ultrasonic  bath  for  about 
half  an  hour.  It  is  then  blow  dried  inside  a  class- 100  clean  hood,  rinsed  with  the  test 
fluid  and  dried  in  an  oven  at  150°C.  The  substrate  is  dipped  into  dilute  Hydrofluoric 
acid  to  remove  any  built  up  oxide,  rinsed  with  DI  water,  test  fluid  and  then  blow 
dried  inside  a  class- 100  hood  using  pure  and  dry  Nitrogen.  One  point  to  note  here  is 
that  Silicon  always  has  a  native  oxide  of  about  30  A  thickness.  Once  the  cell  parts  and 
substrate  are  dry  they  are  assembled  and  mounted  on  the  scanning  stage.  An 
appropriate  amount  of  test  fluid  is  then  introduced  through  the  feed  port  to  form  the 
meniscus.  The  system  is  left  to  equilibrate  with  the  surrounding  for  at  least  half  a 
day  before  taking  any  measurements.  At  this  point  the  Hamaker  constant  can  be 
measured  to  characterize  the  surface. 

The  cell  can  be  evacuated  to  a  high  vacuum  by  bringing  the  cell  pressure  down 
to  50  milli-torr.  Pulling  a  high  vacuum  inside  the  cell  and  then  breaking  the  vacuum 
using  dry  Nitrogen  may  prove  to  be  an  improvement  over  the  above  described 
cleaning  process.  This  possibility  will  be  explored  at  a  later  stage. 

EXPERIMENTAL:  COMPLEMENTARY  CLEANING  STUDIES 

To  investigate  the  effects  of  cleanliness  and  in  the  process  define  the  term  itself 
more  precisely,  we  designed  a  simpler  glass  cell  which  can  be  cleaned  more 
thoroughly  and  which  allows  the  effect  of  different  cleaning  steps  to  be  monitored 
and  controlled.  The  quartz  glass  cell  was  constructed  with  a  ground  glass  joint  at  one 


end,  so  that  a  dessicant  section  can  be  attached  to  preferentially  adsorb  the  water 

vapor.  The  need  for  windows  and  O-rings  was  eliminated  so  that  the  simple  glass  cell 

can  be  easily  cleaned  and  baked.  A  liquid  purification  and  filtration  column  can  also 
be  attached  for  feeding  purified,  filtered  liquid  into  the  cell.  The  liquid  can  be 

purified  by  passing  it  through  a  column  of  silica-gel  and  alumina  to  remove  the 

aqueous  and  polar  impurities  and  then  through  a  membrane  filter  with  a  pore  si/e  of 
0.2  pm. 

Briefly  the  following  cleaning  procedure  will  be  evaluated:  All  glass  parts 
including  the  cell,  dessicant,  feed  and  purification  sections  will  be  rinsed  in  an 
ultrasonic  chamber  with  soap  solution  followed  by  repeated  rinsing  with  purilied 
distilled  water.  They  will  be  dried  in  a  stream  of  pure  filtered  Nitrogen,  rinsed  with 
pure  filtered  fluid,  dried  and  then  baked  at  2 50°C  for  1  hr.  The  wafer  will  be  treated 
with  dilute  HF  (1:100)  and  then  rinsed  with  deionized  water,  dried  with  Nitrogen, 
rinsed  with  pure  filtered  liquid,  dried  and  then  baked  with  the  cell.  The  dessicant 
section  will  be  attached  during  the  baking.  The  cell  will  then  be  cooled  with  the 
dessicant  section  still  attached.  Then  the  cell  will  be  irradiated  bv  short  wavelength 
IJV  light  (184  nm  and  257  inn)  to  clean  the  Silicon  surface  of  adsorbed  organic 
impurities.  This  treatment  is  presumed  to  leave  behind  a  clean,  oxidized  surface. 
Purified,  filtered  liquid  can  now  be  introduced  into  the  cel).  The  importance  and 
effect  of  various  cleaning  variables  will  be  investigated.  Further  cleaning  procedures 
like  creating  an  inert  atmosphere  or  vacuum  in  the  cell  are  also  under  consideration. 
All  these  cleaning  steps  will  be  performed  in  a  controlled  dust  free  atmosphere. 

EmRlMJMS 

In  this  section  we  describe  how  we  measure  the  different  quantities  in  Eq.  (1). 
namely  5(x),  A  experimentally.  The  film  thickness  profiles  of  the  liquid  on  the 
substrate  are  measured  using  interferometry  and/or  ellipsometry.  Both  of  these  two 
methods  are  well  established  (J.G.  Truong  and  P.C.  Wayner,  Jr.,  Cheni.  Rhys.,  &7.  pp. 
4180-4188,  1987),  so  a  detailed  discussion  is  not  included  herein.  We  have 
successfully  tested  the  performance  of  our  cell  with  these  two  optical  techniques.  S(x) 
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i.,  measured  as  a  function  of  power  input  and  fluid  properties.  A  can  he  determined 
by  analyzing  the  ultra-thin  film  section  under  isothermal  condition.  Together  the> 
allow  us  to  calculate  the  integral  heat  sink.  Q.  for  the  region  o(\)  >  0  (from  Fq.  (h).  An 
independent  measurement  of  Q  for  the  total  meniscus  region  is  obtained  by 
measuring  the  total  amount  of  fluid  evaporated  over  a  given  period  of  time.  These 
two  sets  of  data  will  be  analyzed  in  light  of  each  other. 

The  temperature  profiles  will  be  obtained  using  the  thermocouples  attached  to 

A  A  h  '/A  h 

the  hottom  of  the  substrate.  The  interline  heat  flow  parameter.  ,  and  will  also 

v  v 

be  evaluated  and  their  effect  on  fluid  flow  and  heat  transfer  in  the  contact  line  region 
studied.  From  the  film  thickness  profile  we  can  determine  the  evaporative  mass  flux 
leaving  the  surface. 

TUERQRF/OCAL 

One  of  the  major  theoretical  results  to  date  is  Ft].  (37)  in  Appendix  A  entitled 
"The  Effect  of  Interfacial  Mass  Transport  on  Flow  in  Thin  Liquid  Films".  The  equation 
includes  the  combined  effects  of  capillarity  (therefore,  the  capillary  pumping 
parameter),  disjoining  pressure  (therefore,  the  interline  heat  flow  parameter),  gravity 
and  temperature  on  fluid  flow  and  heat  transfer  in  the  contact  line  region.  Since  both 
the  interline  heat  flow  parameter  and  the  capillary  pumping  parameter  are  present  in 
this  equation  their  complementary  effects  can  be  evaluated.  When  the  right  hand 
side  and  the  left  hand  side  of  Fq.  (37)  are  not  equal,  the  terms  in  the  equation  can  be 
used  with  Fq.  (1)  in  the  same  paper  to  analyze  the  transient  case.  Therefore  we  have 

been  able  to  formulate  a  comprehensive  equation  for  general  use.  A  paper  detailing 

numerical  results  based  on  a  simpler  version  of  Fq.  (37)  (i.e..  a  horizontal  film  in 
which  gravitational  effects  are  neglected)  is  being  proposed  for  presentation  at  the 
9th  International  Heat  Transfer  Conference.  Appendix  B  is  a  paper  entitled  "Heat 
Transfer  and  Fluid  Flow  in  an  Evaporating  Extended  Meniscus".  A  note  based  on  onlv 
the  disjoining  pressure  portion  of  this  equation  entitled  "A  Dimensionless  Number  for 
the  Contact  Fine  Ileal  Sink"  was  published  by  the  ASMF  Journal  of  llcat  Transfer  tVol 
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Ill,  813-814,  1989).  Finally,  we  note  that  since  the  theoretical  analysis  includes 
gravity  and  film  tension  (along  with  surface  tension)  the  results  are  very 
comprehensive. 


Further  use  of  this  equation  will,  of  course,  be  directed  at  evaluating  the  results 
of  our  experimental  program. 

CURRENT  STATUS  OF  PRESENTATIONS  AND  PUBLICATIONS 

1  )  Wayner,  P.C.,  Jr.,  '  The  Effect  of  Interfacial  Mass  Transport  on  Flow  in  Thin 
Liquid  Films”,  presented  at  the  International  Symposium  on  Thin  Solid  and  Liquid 
Films.  University  of  Bristol,  England.  11-13  September  1989.  Preliminary  copy 
attached;  to  be  published  in  Colloids  and  Surfaces. 

2  )  Wayner,  Jr.,  P.C.,  Jr.  and  J.  Schonberg,  "Heat  Transfer  and  Fluid  Flow  in  an 
Extended  Meniscus."  to  be  presented  at  the  9th  International  Heat  Transfer 
Conference,  Israel,  August  1990.  Preliminary  copy  attached. 

3)  Wayner,  P.C..  Jr.,  "A  Dimensionless  Number  for  Heat  Transfer  in  the  Contact  Line 
Region."  ASME  ./.  of  Ileal  Transfer ,  Vol.  111.  1089.  PP-  813-814. 

4)  Sujanani.  M.  and  P.C.  Wayner,  Jr.,  "Spreading  at  the  Solid-Liquid-Vapor 
Interline,"  paper  presented  at  the  1  Oth  Symposium  on  Thermophysical  Properties. 
National  Bureau  of  Standards.  Gaithersburg,  MI),  June  20-23,  1988.  (This  presentation 
was  concerned  with  the  characterization  of  the  substrate  surface  and  simple  fluid 
flow  models.) 


8 


QLuniNun  top  °LflTP 


TEST  AREA 


THICKER 
CAPILLARY  GAP 


SILICON  SUBSTRATE 


THINNER  CAPILLARY  GRP 


HEAT  SOURCE 


FIG. - 1  CIRCULAR  CAPILLARY  FEELER 


9 


APPENDIX  A 

THE  EFFECT  OF  INTERFACIAL  MASS  TRANSPORT  ON  FLOW  IN  THIN  LIQUID  MIMS 
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International  Symposium  on  Thin  Solid  and  Liquid  Films 
University  of  Bristol,  Bristol,  England 
11-13  September,  1989 
(Revised  2/1990) 


THE  EFFECT  OF  INTERFACIAL  MASS  TRANSPORT  ON  FLOW  IN  THIN  LIQUID  FILMS 

Peter  C.  Wayner,  Jr. 

Department  of  Chemical  Engineering 
Rensselaer  Polytechnic  Institute 
Troy,  NY  12180-3590,  USA 


The  equilibrium  thermodynamics  and  the  dynamics  of  flow  in  an  extended  meniscus 
that  includes  both  the  relatively  thick  film  region  controlled  by  capillarity  anti  the  ultra- 
thin  film  region  controlled  by  the  liquid-solid-vapor  interfacial  force  field  are  analyzed. 
The  effects  of  capillarity,  disjoining  pressure,  gravity  and  temperature  on  equilibrium, 
fluid  flow  and  mass  transfer  associated  with  the  liquid  film  are  discussed.  The  non¬ 
equilibrium  processes  of  change-of-phase  heat  transfer  and  fluid  flow  in  stationary  and 
in  spreading  thin  films  are  intrinsically  connected  because  of  their  common  dependence 
on  the  intermolecular  force  field  and  gravity.  Criteria  to  determine  the  relative 
importance  of  fluid  flow  and  interfacial  mass  transfer  based  on  the  film  profile  and  the 
thermophysical  properties  of  the  system  are  developed.  An  evaporation/condensation 
mechanism  is  found  to  be  important  in  contact  line  motion. 
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dlnf 


dyf 

RTSAn 


As 

RTAn 


dT 


(9) 


Derjaguin  and  Zorin  |5]  essentially  used  Eq.  (9)  to  study  adsorption  on  a  superheated 
solid  surface  placed  close  to  but  slightly  above  a  pool  of  liquid.  Using  Eq.  (9)  with  dlnf  = 
0  and  Eq.  (7)  for  dyf  gives 

dn  =  (sv  -  s/)dT  (10) 


Integration  of  Eq.  (10)  over  a  small  temperature  change  (Tv,  n  =  0;  T/v,  n)  during  which 
(sv  -  s /)  is  approximately  constant  leads  to 

n  =  (sv  -  s /)  (T/v  -  Tv)  (11) 


Using  n  =  -A/6tr53  for  film  thicknesses  less  than  20  nm  and  (sv  -  s /)  =  Ah/Tj,  where  T,  = 
(T/v  +  Tv)/2, 


-  A  T  i  i  ,  t 

§  =  ( -  )l/3 

WT/v  -  Tv)Ah; 


(12) 


Eq.  (12)  can  be  used  to  calculate  the  thickness  of  an  adsorbed  layer  of  liquid  on  a 

superheated,  (T/v  -  Tv),  solid  surface  from  the  value  of  the  heat  of  vaporization  of  the 

liquid  film,  Ah,  and  the  Hamaker  constant,  A.  Conversely,  the  Hamaker  constant  cam  be 
obtained  by  measuring,  8(T/V,  Tv). 

For  small  changes  in  the  fugacity  in  a  range  where  the  fugacity  can  be  replaced  by 

the  vapor  pressure,  Eqs.  (7  and  9)  give  on  integration  with 

(An)'1  =  -V/ 


V/nPv  PvV/Ah 

P/V  -  Pv  =  -  RT/v  +  RT/vTv  (T/v  '  Tv)  (13) 

In  this  equation,  Pv  is  a  reference  vapor  pressure  of  a  thick,  flat,  liquid  film, 

II  — >  0,  with  a  surface  temperature  T/v  =  Tv  and  Ah  is  the  ideal  heat  of  vaporization  of  the 
film.  Equation  (13)  can  be  modified  to  include  curvature  at  the  liquid  vapor  interface  by 

replacing  n  by  n  +  o/vK  in  which  K  is  the  curvature  and  o/v  is  the  surface  tension  of  the 
liquid-vapor  interface: 

V/Pv  PyV/Ah 

P/V  -  Pv  =  -  (n  +  c/v K)  +  R~T~fv~  (T/v  -  Tv)  (14) 

Equation  (14)  also  represents  the  change  in  equilibrium  vapor  pressure  relative  to  the 
pool  surface  of  the  extended  meniscus  presented  in  Figure  (1)  for  a  change  in 
temperature  of  (T/v  -  Tv)  and  hydrostatic  pressure  equal  to  p^gx  =  (fl  +  o/vK).  This  is 
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approximate  since  buoyancy  has  been  neglected.  An  integration  in  the  vapor  space  ai 
constant  Tv  gives  for  small  changes  in  the  vapor  pressure 


P 


VX 


-  Pv 


MgPyX 

RTV 


(15) 


Combining  Eqs.  (14  and  15)  with  Pvx  =  P/v  gives  for  the  equilibrium  case 
V/Pv  _  MgPv  PvV/Ah 

P/v  -  Pvx  =  -  (n  +  o/vK)  +  v  x  +  Rf (T/v  -  Tv)  =  0  (16) 

lor  T/v  =  Tv,  this  reduces  to  the  identity 

V/(H  +  o/vK)  =  V^gx  =  (Mgx)/  =  (Mgx)v  (17) 


If  Pvx  *  P/v,  Eq.  (16)  can  be  used  along  with  kinetic  theory  to  calculate  the  rate  of 
evaporation  from  (or  condensation  on)  a  curved  thin  film  as  demonstrated  in  the  next 
section. 

Interfacial  mass  flux. 


In  1953  Schrage  [  1 7 1  reviewed  the  literature,  presented  and  discussed  the 
following  equation  based  on  kinetic  theory  relating  the  net  mass  flux  of  matter  crossing 
a  liquid-vapor  interface  to  a  jump  change  in  interfacial  conditions  at  the  interface: 


Herein,  we  presume  that  the  net  mass  flux  crossing  the  interface  (e.g.,  evaporation) 
results  from  a  small  vapor  pressure  drop  across  an  imaginary  plane  at  the  interface  in 
which  P/v,  is  the  quasi-equilibrium  vapor  pressure  of  the  liquid  film  at  (T/v,  K,  H.  \)  and 
Pvx  is  the  equilibrium  vapor  pressure  of  a  reference  bulk  liquid  (K  =  0,  FI  =  0,  x)at  a 
temperature  Tv.  Neglecting  resistances  in  the  bulk  vapor  space,  Pvx  and  Tv  can  exist  at  a 
short  distance  from  the  interface,  and  a  resistance  to  evaporation  at  the  interface  can  be 
defined  using  Eq.  (18).  Using  T1^2  =  T1^2,  this  can  be  rewritten  as 


M 

m  =  Cl  <2«RTi)l/2  (P,V  '  P'"‘> 


(19) 


Wayner,  et  al.  Ill]  used  an  extended  Clapeyron  equation  for  the  variation  of  equilibrium 
vapor  pressure  with  temperature  and  disjoining  pressure  in  a  horizontal  thin  film  to  obta 
Eq.  (20)  for  the  vapor  pressure  difference  in  Eq.  (19). 


P/v 


VO 


P  vM  Ahm 

RTvT/v 


(T/v  -  Tv)  + 


V/Py 

RT/v 


(P/  -  Pv) 


(20) 
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We  note  that,  when  interfacial  effects  are  important,  the  effective  pressure  in  the  liquid, 
P/,  is  not  necessarily  equal  to  the  pressure  in  the  vapor  or  to  its  normal  vapor  pressure. 
For  example,  the  equilibrium  vapor  pressure  of  a  small  droplet  is  not  the  same  as  that  ol 
a  flat  pool  of  bulk  liquid  at  the  same  temperature,  and  the  pressure  inside  the  droplet  is 
greater  than  the  surrounding  pressure  because  of  surface  tension.  The  concept  of 
pressure  is  further  clouded  by  the  concept  of  internal  pressure  due  to  cohesion  which  is 
extended  in  thin  films  to  include  adhesion.  Therefore,  P/v  is  used  to  designate  the 
equilibrium  vapor  pressure  of  an  interface  at  T/v  changed  by  interfacial  forces  and  (P /  - 
Pv)  designates  the  pressure  change  in  the  liquid  due  to  interfacial  forces.  For  example, 
the  above  equation  can  be  used  to  calculate  the  required  subcooling  needed  to  achieve 
equilibrium  between  a  curved  interface  (droplet)  and  a  flat  interface  when  the  vapor 
pressure  is  changed  by  surface  tension. 

Combining  Eqs.  (19  &  16)  gives  for  the  non-equilibrium  case 


m  =  Ci  ( 


M 


2ttRTi 


)  1/2 


P  V  M  A  h  r 
RTVT/V 


V/Pv  MgPv 

(t,v ' Tv) '  rt/7  (n  +  °/vK)  +  “rt7  x)} 


(21) 


which  can  be  used  to  calculate  the  interfacial  mass  flux. 


Herein,  there  are  two  interfacial  effects  that  can  cause  an  effective  pressure  jump  tit 
the  liquid-vapor  interface:  capillarity,  o/vK,  and  disjoining  pressure.  11.  The  significance 
of  the  disjoining  pressure  herein  is  that  the  vapor  pressure  of  an  adsorbed  completely 
wetting  liquid  film  is  reduced  by  interfacial  forces  and  therefore  a  superheated  adsorbed 
liquid  film  can  exist  in  (  vapor  pressure)  equilibrium  with  a  bulk  liquid  at  a  lower 
temperature.  For  convenience,  we  will  emphasize  spreading  (zero  contact  angle)  fluids 
and  will  use  the  following  sign  convention 

P/  -  Pv  =  -o/vK  -  fl  (22) 


_B  _  A 
§n  “  '  83 


(23) 


The  approximation,  n  -  3,  in  Eq.  (23)  restricts  its  use  to  thicknesses  5  ^  20  nm.  However, 
in  the  numerical  examples  presented  below  this  is  acceptable  because  either  5  jf  20  nm 
or  fl  r  0.  The  equations  can  be  easily  modified  for  general  use  and/or  for  other  cases.  It 
is  useful  to  make  the  mass  flux  dimensionless  using  an  ideal  mass  flux  based  on  the 
temperature  change  only: 


m 


id  = 


Ci  ( 


M 


2rcRT 


)l/2 


P  v  M  A  h  m 

(T777)  (T'V ' Tv) 


(24) 


M  = 


m 

m icl 


Mi 

M AhmAT 


(o/v K  +  n  -  p/gx) 


(25) 
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with  AT  =  (T/v  -  Tv).  Taking  IV1  =  0,  n  =  3,  B  =  A  =  A/6 k,  K  =  0  in  Eq.  (25),  Eq.  (12)  for  the 
thickness  of  a  flat  adsorbed  liquid  film,  80,  on  a  superheated  surface  with  AT0  =  (T/v  -  Tv) 
can  be  recovered. 


Restricting  the  following  material  to  a  constant  interfacial  temperature  difference,  ATU  = 
('i'/v  -  Tv),  and  defining  a  reference  disjoining  pressure,  n0  3  -  A/S^,  Eq.  (25)  becomes 


m  =  i  -  (o/vk  +  n  -  p/gx)  rro' 


(26) 


A  constant  temperature  difference  simplifies  the  analysis  and  should  not  alter  the 
significance  of  the  results.  It  is  useful  to  keep  in  mind  that  n0  could  be  replaced  by  the 
equivalent  temperature  difference  using  Eq.  (12). 


Using  rj  =  8/80  and  E,  =  x/x0  the  curvature,  K,  for  small  interfacial  slopes,  can  be  written  as 
8" 


_ _ _  n  8o  d2r| 

K  =  [1  +  (8')2]3/2  "  5  =  x2  d^2 


(27) 


Defining  the  parameter  3  Nc  3  Oyv80/n0x^  with  n0  =  r|3n  and  p/gx  =  n0X  gives 


M  =  1  -  3Ncti"  -tt3  +  X  ( 2 X ) 

With  dimensionless  terms  the  prime  refers  to  differentiation  with  respect  to  E.  otherwise 
\.  Eq.  (28)  gives  the  change  in  the  dimensionless  evaporation  rate  resulting  from  the 

effects  of  gravity  and  interfacial  forces  on  the  vapor  pressure.  Using  aAT0  -  mAhm,  this 
can  also  be  written  as  a  dimensionless  liquid-vapor  interfacial  heat  transfer  coefficient: 

=  I  -  3Ncr|"  -  q-3  +  X  (29) 

aul 

Therefore,  an  interfacial  heat  transfer  coefficient  can  be  obtained  from  first  principles. 

In  1 11]  only  the  changes  due  to  disjoining  pressure  and  temperature  were  analyzed. 
Considerable  insight  can  be  obtained  from  Eq.  (28).  We  find  that  at  the  leading  edge  of 
an  evaporating  horizontal  (X  =  0,  due  to  g  =  0)  flat  (rj"  =  0)  thin  film,  rj  =  1,  the 
evaporation  rate  can  be  zero  even  through  the  superheat  can  be  substantial.  This 
condition  is  required  to  define  a  stationary  interline  for  an  evaporating  thin  liquid  film 
with  varying  thickness  which  has  an  equilibrium  contact  angle  equal  to  zero.  The 
interline  has  been  defined  as  the  location  of  the  junction  of  an  evaporating  thin  film  and 
a  non-evaporating  adsorbed  film  (11].  However,  when  the  dimensionless  film  thickness 
is  11  =  4.6  (X  =  0,  tj "  =  0),  the  evaporation  rate  will  be  99%  of  its  substantial  ideal  value. 
When  rj  <  1 ,  condensation  occurs.  These  equations  allow  the  interfacial  forces  and 
thermal  effects  to  be  compared  and  combined,  and  have  many  uses  in  evaporating. 
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condensing  and/or  equilibrium  systems..  For  curvature  to  have  the  same  effect,  use  of 
Eqs.  (12  and  26)  with  n  =  0  and  g  =  0  gives 


K0 


AhAT0 

°/vTv 


(30) 


Using  F.q.  (12)  and  n0  =  Pygx,  the  following  equivalents  are  obtained  for  octane  on  SiO  ->  at 
208. 16K  with  A/6tt  =  -3.18  x  10-22  j. 


A  l'o,  K  K0,  m1  n0,  N/m2  80,  nm  x,  m 

0.1  3.94  x  106  8.49  x  1Q4  1.55  12.4 


This  demonstrates  that  the  decrease  in  film  thickness  associated  with  a  A  T  =  0. 1  K  is 
equivalent  to  the  pressure  difference  associated  with  a  static  liquid  head  of  12.4  m. 

This  is  a  pressure  drop  available  for  fluid  flow. 

Fluid  Mechanics. 

Fluid  flow  in  a  curved  thin  film  is  also  controlled  by  the  "pressure"  gradient.  Using 
the  definitions  for  q,S,,  x0,  n0  and  the  small  slope  assumption  given  above  leads  to 

d  3n„  tV  X’ 

dl (°/vK  +  n "  P^x)  =  77  v  *  7>  <3  1  > 

Using  the  lubrication  approximation  while  neglecting  surface  shear  stress,  the  mass  flow 
rate  per  unit  width  in  a  liquid  film  can  be  approximated  by 

r  =  7  dTP/  +  P/S*l  (32) 


Using  Eqs.  (22,  31  and  32)  while  assuming  pv  «  p /  gives 


r  = 


§3n0 

x0v 


(Net) 


tV 

Tl4 


(33) 


and  the  following  equation  for  the  dimensionless  mass  flow  rates  in  the  film 


r  s 


rx0v 

8ono 


=  Ncrj^r)"' 


tV 

r\ 


n-3x' 

3 


(34) 


DISCUSSION 


Using  Eq.  (35)  for  a  steady  state,  stationary  thin  film  with  phase  change 


dr 


(35) 
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and  Eq.  (36)  with  Eq.  (28)  gives  Eq.  (37) 


dT  m_xjv 
=  A 


in  which  x^  =  -A/v  m  1(1 
n'  r)3x' 

(—  +  ~~  -  Net)  V’V  =  l  -  3  Net) ”  -  r|-3  +  X 
t|  3 


(36) 


(37) 


Using  four  bondary  conditions,  Eq.  (37)  can  be  numerically  solved  to  determine  the  heat 
transfer  characteristics  of  a  stationary,  steady  state,  thin  film.  When  the  right-hand  side 
and  the  left-hand  side  of  Eq.  (37)  are  not  equal,  the  above  material  can  be  used  with  Eq. 
(1)  to  analyze  the  transient  case. 

I 

An  evaluation  of  Eq.  (21)  demonstrates  that,  in  an  isothermal  horizontal  ultra-thin 
film,  evaporation  or  condensation  occurs  unless  (n  +  o/vK)  =  0.  Using  simple  models,  this 
case  without  phase  change  has  been  previously  discussed  for  the  leading  edge  of  a  film 
(Tl  >  0,  K  <  0)  by  deGennes  and  Joanny  [18]  and  for  the  partial  wetting  case  (n  <  0,  K  >  0) 
by  Wayner  [19].  The  results  demonstrate  that  the  interface  becomes  planar  very 
ra pidly. 

For  motion  to  occur  in  an  isothermal,  ultra-thin,  horizontal  film,  Eq.  (32) 
demonstrates  that  (n  +  o/vK)  *  0.  This  indicates  that  fluid  flow  due  to  interfacial  forces 
in  an  "isothermal"  thin  film  must  be  associated  with  change-of-phase  heat  transfer  (Eq. 
(21)).  Therefore,  at  the  leading  edge  of  a  spreading  ultra-thin  film  a  small  temperature 
gradient  due  to  change-of-phase  heat  transfer  is  present.  However,  this  would  be 
extremely  small  with  low  volatile  fluids.  In  general  the  temperature  gradient  would  not 
effect  the  above  analysis  unless  a  relatively  large  gradient  was  applied. 

RELATIVE  IMPORTANCE  OF  INTERFACIAL  MASS  FLUX  TO  LIQUID  FLOW  RATE  GRADIENT 

Draining  Vertical  Thin  Film 

Neglecting  curvature  and  curvature  gradient  terms  as  compared  to  gravity 
for  the  case  of  a  relatively  thick  (n  =  0),  isothermal  (AT  =  0)  draining  film,  the  ratio  of  Eq. 
(21)  to  the  derivative  of  Eq.  (32)  leads  to 


m  _  -EV/PyVx 
T  ’  =  RTV828‘ 

in  which 


E  =  Ci  ( 


M 

2nRTj 


•)l/2 


(38) 


(39) 
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For  8  =  1(F6  m,  x  =  1(F2  m  with  octane  at  298  K  (Pv  =  14  mmHg,  Ev  =  4  x  1(F9  m).  we  Find 


m_  _  4.8(10:3) 
r  '  S' 

Recent  studies  of  draining  liquid  films  of  octane  found  that  8'  =  6  x  1(F3  at  x  =  10  2  m 
when  8  =  10-6  m  [20].  Therefore,  both  evaporation  and  liquid  flow  would  effect  the  rate 
of  decrease  of  film  thickness  for  these  conditions.  In  these  studies,  the  films  were 
observed  to  drain  faster  than  that  anticipated  due  to  gravitational  forces  only. 

Additional  interfacial  effects  due  to  (o/vK)'  were  also  probably  present.  Further,  this 
ratio  and  those  discussed  below  do  not  include  the  resistance  to  diffusion  in  the  vapor 
space. 

Spreading  Horizontal  Film  in  the  Region  5  -»  0  (Contact  Line  Motion) 

For  the  case  of  a  horizontal  (g=0)  isothermal  (AT  =  0),  spreading  (H  =  -A/83), 
slightly  tapered  ultra-thin  film,  Eqs.  (21,  22,  and  32)  with  n  »  aK  and  (oK)’  =  0  gives 


m  _  EV/Pvv 

I'  1  ~  RT/V8[58"  -  (61)2! 


(4  1  ) 


For  octane  at  298  K  with  8  =  10'9  m,  8’  -10*6,  8"  ~10‘3  nr1,  [88"  -  (8')2]  is  of  the  order  1(F 

'2,  and 


m 

r 


4.9  x  10-13 
(10-9)  (10-12) 


~  0(108) 


(42) 


It  appears  that  contact  line  motion  is  primarily  due  to  an  evaporation/ 
condensation  mechanism  for  these  conditions.  Once  a  change  of  phase  occurs, 

AT  t  0,  and  the  value  of  the  ratio  in  Eq.  (42)  would  decrease.  However,  this  would  not 
affect  the  basic  conclusion.  Although  an  additional  resistance  due  to  diffusion  would  be 
present  with  a  non-condensing  second  component  in  the  vapor,  the  diffusion  path  length 
can  be  extremely  small.  We  note  that  Hardy  and  Doubleday  found  that  primary  film 
spreading  was  a  function  of  molecular  weight  and  vapor  pressure  when  the  fluids  had  a 
sensible  vapor  presure  |2|. 

Contact  Line  Heat  Sink 

Using  T  =  pSV  to  define  an  average  velocity  in  the  film,  Eq.  (32)  and  the 
above  assumptions  (g  =  0,  K  =  0)  lead  to 


V  = 


A8' 

p82 


(43) 
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Using  the  value  8'  =  -10'6  obtained  above  with  A  =  -3.18  x  1022  J  (octane  -  SiCb),  p  - 
5.16  x  10‘4  N-s/m2  and  8  =  10‘9  m  in  Eq.  (43)  gives  V  =  6.1  x  10'7  m/s  for  perspective. 

Assuming  that  all  the  liquid  which  flows  through  a  plane  perpendicular  to  the 
substrate  at  a  given  location  evaporates  in  the  thinner  portion  of  a  steady  state, 
stationary  film  with  AT  >  0,  a  contact  line  heat  sink  can  be  defined  as 


Q  =  TAhm 


A8  Ahm 
v8 


(44) 


Therefore,  it  is  possible  to  obtain  the  contact  line  heat  sink  knowing  the  thermophysical 
properties,  the  maximum  thickness  of  the  region  and  the  slope  at  this  thickness.  The 

slope  is  a  measure  of  the  departure  from  equilibrium  leading  to  fluid  flow;  the 
temperature  difference  and  reference  condition  indicate  whether  or  not  phase  change 
occurs. 

Spreading  Horizontal  Film  (Incomplete  Wetting) 

The  following  simple  model  has  been  used  to  describe  the  contact  line  region 
of  an  equilibrium  liquid  film  with  a  finite  contact  angle  [19], 

n  +  a/vK  =  0  (45) 

In  this  case  n  <  0  and  K  >  0,  or  we  might  say  that  the  disjoining  pressure  increases  the 
vapor  pressure  and  the  curvature  decreases  the  vapor  pressure.  Since  fluid  motion  in 
the  liquid  film  would  result  from  an  increase  in  the  curvature  gradient,  this  curvatme 
change  would  also  lead  to  condensation  and  a  change  in  the  apparent  contact  angle.  1  he 
processes  can  be  separated  only  by  degree  of  importance. 

CONCLUSION 

The  physical  phenomena  associated  with  spreading  have  been  extensively 
addressed  using  various  models  for  low  vapor  pressure  fluids  (e.g.,  see  the  reviews  by 
deGennes  [21],  Cazabat  [22],  and  Neogi  and  Miller  [9]).  They  range  from  the  rolling 
motion  observed  by  Dussan  and  Davis  [23]  to  the  surface  diffusion  model  discussed  by 
Ruckenstein  and  Dunn  [24],  As  demonstrated  above,  change  of  phase  heat  transfer  is 
also  important  for  many  conditions.  Bascom,  et  al.  [6]  citing  deBoer  [25 1  noted  that 
surface  diffusion  can  be  energetically  more  favorable  than  evaporation.  Indeed,  tie  Boer 
125 1  discussed  this  situation  for  the  case  of  a  two-dimensional  gas.  He  also  noted  that 
lateral  migration  would  be  restricted  when  two-dimensional  condensation  occurred. 
Additional  research  on  spreading  as  a  function  of  film  thickness  is  warranted.  In 
general,  the  non-equilibrium  processes  of  fluid  flow  and  change  of  phase  heat  transfer 
are  intrinsically  connected  because  of  their  common  dependence  on  the  intermolecular 
force  field  and  gravity.  The  relative  importance  of  the  various  mechanisms  depend  on 
the  physical  properties  of  ultra-thin  films  and  how  they  are  effected  by  pressure,  ( i . c . , 
cohesion,  adhesion,  thickness,  curvature)  and  temperature. 
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NOMENCLATURE 


A  =  Hamaker  constant  (6tcA) 

C  =  condensation  coefficient 

f  =  fugacity 

g  =  acceleration  due  to  gravity 

h  =  enthalpy/  volume,  heat  transfer  coefficient 

K  =  curvature 

M  =  molecular  weight 

in  -  interfacial  mass  flux 

n  =  molar  density 

P  =  pressure 

R  =  gas  constant 

s  =  entropy/volume 

T  =  temperature 

t  =  time 

V  =  molar  volume 

X  =  dimensionless  position 

a  =  interfacial  heat  transfer  coefficient 

A  =  difference 

o  =  film  thickness 

I  =  mass  flow  rate  /  width 

y  =  film  tension 

il  =  dimensionless  thickness 

li  =  chemical  potential/mole,  dynamic  viscosity 

v  =  kinematic  viscosity 

II  =  disjoining  pressure 

p  =  density 

a  =  interfacial  tension 

q  =  dimensionless  position 

stt  bscripts 

g  =  includes  gravity 

i  =  average  value 

I  =  liquid 

lv  =  liquid-vapor  interface 

m  =  unit  mass 

T  =  tangential 

v  =  vapor 

x  =  at  x 

superscripts 

ill  =  ideal 

'  ,  "  =  derivatives 
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APPENDIX  B 

HEAT  TRANSFER  AND  FLUID  FLOW  IN  AN 
EVAPORATING  EXTENDED  MENISCUS 
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HEAT  TRANSFER  AND  FLUID  FLOW 
IN  AN  EVAPORATING  EXTENDED  MENISCUS 


P.C.Wayner.Jr.  and  J.  Schonberg 


Department  of  Chemical  Knginecring 
Rensselaer  Polytechnic  Institute 
Troy.  N  Y.  121X0 


The  rate  of  evaporation  from  a  steady  extended 
meniscus  formed  between  the  exits  of  two  facing  slot 
feeders  is  analyzed.  The  meniscus  is  symmetric  about 
the  ccnterplane.  The  physicochemical  phenomena 
associated  with  cohesion  and  adhesion  at  interfaces 
arc  used  to  formulate  transport  equations  unique  to 
evaporating  thin  films  which  arc  then  numerically 
solved  for  two  types  of  menisci.  One  type  involves  an 
adsorbed  thin  film  dominated  by  disjoining  pressure 
effects,  a  thicker  portion  dominated  by  capillary 
effects  and  a  transition  region  between.  The  second 

type  has  no  thin  film  portion:  the  transition  region 

straddles  the  plane  of  symmetry.  In  a  sample 

calculation  the  heat  flux  is  found  to  be  roughly  nine 
times  higher  for  the  latter  type  of  meniscus.  For  the 
conditions  studied,  we  find  that  a  smaller  value  of  the 
llamaker  constant  (weaker  adsorption)  leads  to  an 
Increase  in  the  relative  importance  of  capillarity  and 
to  a  higher  evaporation  rate. 


I  INTRODUCTION 

Fluid  flow  and  evaporation  in  an  extended 
meniscus  that  includes  both  a  relatively  thick  region 
controlled  by  capillarity  and  an  ultra-thin  region 
controlled  by  the  liquid-solid-vapor  interfacial  force 
field  (disjoining  pressure)  is  of  importance  to  many 
heat  transfer  processes.  When  the  film  thickness  is 

less  than  10’7  m,  long  range  intermolccular  forces 
start  to  affect  the  pressure  distribution  in  the  liquid 
and  its  vapor  pressure.  For  fluids  that  completely  wet 
the  substrate,  these  effects  are,  of  course,  extremely 

.7 

small  for  thicknesses  of  the  order  of  10  m.  However, 
they  become  the  predominant  effect  as  the  film 
becomes  ultra-thin.  Herein,  we  evaluate  these 
complementary  effects. 

Dcrjaguin  cl  al.  (c.g.  1957,  1965)  laid  the 
foundation  for  the  study  of  the  mechanics  and 
thermodynamics  of  ultra-thin  films  when  they 
measured  the  thickness  of  an  adsorbed  liquid  film  on 
a  superheated  substrate,  discussed  the  flow  of  an 
evaporating  ultra-thin  film  in  a  capillary  and  defined 
the  disjoining  pressure  to  account  for  the  modified 
force  field  in  the  liquid  Potash  and  Wayncr  (1972) 
used  these  concepts  to  describe  evaporation  from  a 
stationary,  steady-state,  evaporating  extended 
meniscus  on  a  vertical  flat  plate.  Miller  (1973)  studied 


the  stability  of  moving  liquid-vapor  interfaces  as  a 
result  of  phase  change.  Wayncr  ct  al.  (1976)  used  an 
extended  Clausius-Clapcyron  equation  to  define  an 
interline  heat  transfer  coefficient.  In  1979  these 
concepts  were  used  to  model  evaporation  from 
grooves  (Holm  and  Goplen  (1979))  .  Moosman  and 
Homsy  (1980)  used  perturbation  theory  to  describe  the 
profile  change  of  an  evaporating  extended  meniscus 
relative  to  the  static  isothermal  profile.  The  thermal 
resistance  of  the  liquid  itself  was  included.  The  heat 
flux  was  found  to  be  a  maximum  in  the  transition 
region.  In  the  capillary  region  (for  the  particular 
parameters  chosen)  the  film  resistance  hindered 
evaporation.  In  the  thin  adsorbed  film  disjoining 
pressure  hindered  evaporation.  Herein.  film 
resistance  is  neglected;  however  the  fully  non-linear 
equations  arc  numerically  analyzed;  thus  the  heal 
flux  may  be  larger  than  in  the  linearized  process 
(Moosman  and  Homsy.  19X0).  Furthermore  in  this 
work  the  meniscus  need  not  terminate  in  an  adsorbed 
thin  film. 

Studies  by  Parks  and  Wayncr  (1987)  and 
Mirzamoghadam  and  Catton  (1988)  extended  the 
literature  to  include  other  effects  (e.g.  surface  shear) 
Mirzamoghadam  and  Canon  (1988)  considered  an 
evaporating  meniscus  attached  to  a  heated  plate 
partially  immersed  in  a  liquid  pool  at  various  angles 
An  approximate  integral  analysis  was  employed  to 
study  the  whole  process  from  the  adsorbed  thin  film 
to  the  region  of  natural  convection.  Furthermore  the 
influence  of  temperature  on  surface  tension  was 
considered.  Like  Moosman  and  Homsy  (1980)  the 
resistance  of  the  meniscus  to  thermal  conduction  was 
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included.  Das  and  Gaddis  (1987)  evaluated  the 
pressure  distribution  near  an  evaporating  contact 
line.  Herein,  we  apply  these  concepts  to  a  new 
meniscus  configuration  (Fig.  1)  which  naturally  leads 
to  the  evaluation  of  important  new  boundary 
conditions  in  the  numerical  solution  of  the  full  non¬ 
linear  equation  for  an  evaporating  single  component 
film  with  a  constant  surface  temperature..  The  basic 
equations  are  reformulated  so  that  a  dimensionless 
group,  Nc,  represents  the  relative  importance  of 
capillarity  to  a  reference  disjoining  pressure  which, 
in  turn,  is  a  function  of  the  substrate  superheat 


Fquation  (3)  can  be  used  to  calculate  the  required 
subcooling  needed  to  achieve  equilibrium  between  a 
curved  interface  (droplet)  and  a  flat  interface  when 
the  vapor  pressure  is  changed  by  surface  tension 
Combining  Eqs.(2)  &  (3)  gives  the  interfacial  mass 
flux  as 


c  1/2 

M  (2irRT' 


PyMAh 

RTvT/v 


(T/v  -  Tv) 


IP I  -  Pv>} 


(4) 


2  THEORETICAL  BACKGROUND 

2.1  Interfacial  mass  flux:  Schragc  (1953)  reviewed 

the  literature,  and  presented  and  discussed  Eq.  (1) 
based  on  kinetic  theory  relating  the  net  mass  flux  of 
matter  crossing  a  liquid-vapor  interface  to  a  jump 
change  in  intcrfacial  conditions  at  the  interface: 


The  net  mass  flux  across  the  interface  results  from  a 
small  vapor  pressure  drop  across  an  imaginary  plane 
at  the  interface  in  which  P/v  is  the  equilibrium  vapor 
pressure  of  the  bulk  liquid  at  T/v.  and  Pv  is  the 
equilibrium  vapor  pressure  of  the  reference  bulk 
liquid  at  a  temperature  Tv.  Neglecting  resistances  in 
the  bulk  vapor  space.  Pv  and  Tv  can  exist  at  a  short 
distance  from  the  interface,  and  a  resistance  to 
evaporation  at  the  interface  can  be  defined  using  Eq. 
U).  Using  t'^  -  T*^“ .  this  can  be  rewritten  as 


Herein,  there  are  two  intcrfacial  effects  that  can 
cause  a  pressure  jump  at  the  liquid-vapor  interface: 
capillarity.  yK,  and  disjoining  pressure,  n  A  recent 
review  article  (Ivanov,  1988)  describing  the  concepts 
is  available.  The  significance  of  the  disjoining 
pressure  herein  is  that  the  vapor  pressure  of  an 
adsorbed  liquid  film  is  reduced  (complete  wetting 
case)  by  intcrfacial  forces  and  therefore  a 
superheated  adsorbed  liquid  film  can  exist  in  (vapor 
pressure)  equilibrium  with  a  bulk  liquid  at  a  lower 
temperature.  For  convenience,  we  will  only  look  at 
spreading  (zero  contact  angle)  fluids  herein  and  will 
use  the  sign  convention 


P/  -  Pv  =  -yK  -  fl 


(5) 


(6) 


in  which  A  is  the  Hamakcr  constant.  K  the  curvature, 
y  the  surface  tension  and  8  the  film  thickness. 


m  =  C|  <2itRT,l/2  <P/v  '  Pv  1  (2) 

Wavncr  et  al  (1976)  used  an  extended  Clapeyron 
equation  for  the  variation  of  equilibrium  vapor 
pressure  with  temperature  and  pressure  to  obtain  Eq. 
(3)  for  the  vapor  pressure  difference  in  which  Pv/v 
replaces  P/v 


The  approximation,  n  =  3.  used  in  Eq.  (6)  restricts  the 
use  of  Eq.  (6)  to  thicknesses  8  *  20  nm  (Truong  and 

Wayncr,  1988).  However,  the  numerical  examples 
presented  below  demonstrate  that  this  is  acceptable 
for  some  cases  in  which  yK  becomes  dominant  before 
this  thickness.  Further,  the  equations  can  be 
modified  for  general  use  for  other  cases.  It  is  useful 
to  make  the  mass  flux  dimensionless  using  an  ideal 
mass  flux  based  on  the  temperature  change  only: 


,1/2  PyyMjih 


'RTvT/v1  (T,v 


Tv  I  (7) 


in  which  Ah  is  the  heat  of  vaporization  of  the  film 
and  V/  is  the  molar  volume. 


M  —  —  =  1 
m  ’  d 


MAhAT  (yK  +  n) 


Wc  note  that,  when  intcrfacial  effects  arc  important, 
the  "pressure".  P /,  in  the  liquid  is  not  necessarily 
equal  to  the  pressure  in  the  vapor  or  to  its  normal 
vapor  pressure.  For  example,  the  equilibrium  vapor 
pressure  of  a  small  droplet  is  not  the  same  as  that  of  a 
flat  pool  of  bulk  liquid  at  the  same  temperature  and 
the  pressure  inside  the  droplet  is  greater  than  the 
surrounding  pressure  because  of  surface  tension  The 
concept  of  pressure  is  further  clouded  by  the  concept 
of  internal  pressure  due  to  cohesion  which  is 
extended  to  include  adhesion  herein.  Therefore,  Pv/v 
is  used  to  designate  the  equilibrium  vapor  pressure  of 
an  interface  at  T/v  affected  by  intcrfacial  forces. 


with  AT0  =  (T/v  -  Tv).  Taking  M  =  0,  n  =  3.  B  =  A.  K  =  0  in 
Eq.  (8).  the  thickness  80  of  a  flat  adsorbed  liquid  film 
on  a  superheated  surface  can  be  obtained: 


6o 


-V/TvA 

<MAhAT„) 


1/3 


(9) 


Restricting  the  following  material  to  a  constant 
intcrfacial  temperature  difference.  AT0  and  defining 
a  reference  disjoining  pressure.  Flo  * 

-  A/82,  Eq.  (8)  becomes 
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M  =  1  -  (yK  +  n)  n'J 


(10) 


(18) 


Using  i)  =  S/So  xid  4  =  */xo  O’®  curvature.  K,  for  small 
interfacial  slopes,  can  be  written  as 


K  =(1  +  (S')3 ,3/2 


s„  _  s^ 

-  .J  *2 


(ID 


Defining  the  parameter  3  Nc  =  YS0/n0x^  with  n0  =  r|-^n 
gives 


M  =  1  -  3Ncq "  -  rr3  (12) 

Eq.  (12)  gives  the  reduction  in  the  evaporation  rate 
resulting  from  the  effect  of  interfacial  forces  on  the 

vapor  pressure.  Using  aAT0  =  mAh.  this  can  also  be 
written  as  a  dimensionless  liquid-vapor  interfacial 
heat  transfer  coefficient: 

=  1  -  3Nct|"  -tt3  (13) 

a10 


r 


rx0v 

5ono 


=  Ncq3q"  - 


n! 

n 


Using 


dr 

T~  =  -  m 
dx 


and 


dT 

<14  = 


A 


-M 


(19) 


CO) 


in  which  =  -A/v  mia,  with  Eq.  (12)  gives 

-  Nct|3t)'")'  =  I  -  3Ncq”  -  q ‘ 3  (21) 

The  lefthand  side  of  Eq.  (21)  represents  fluid  flow  and 
the  righthand  side  represents  intcrfacial  mass  How. 


Considerable  insight  can  be  obtained  from  Eq.  (12). 
Consider  two  hypothetical  meniscii.  Each  is  a 
perfectly  flat  film,  one  is  deeper  (q  =  4.6)  and  one  is 
shallower  (q  =  1.0).  Each  is  superheated;  however, 
there  is  no  evaporation  from  the  latter  meniscus  due 
to  adsorption.  However  the  evaporation  rate  from  the 
former  meniscus  is  99%  of  its  ideal  value.  For 
curvature  alone  to  prevent  evaporation  it  must  be 
equivalent : 


3  RESULTS  AND  DISCUSSION 

3.1  Numerical:  In  principle  the  descriptive  Eq 

(21)  may  be  solved  as  written,  however,  the  parameter 
Nc  may  be  small.  In  that  case  the  highest  order 
derivative  is  multiplied  by  a  small  parameter  which 
creates  numerical  difficulties.  This  problem  was 
avoided  by  rescaling  the  dimensionless  distance.  A 
new  distance  X  is  thus  defined  according  to 


„  _  MAhATp 

0  ”  YV/Tv 


(  14) 


(22) 


This  follows  from  Eqs.  (9)  and  (10)  and  the 
requirement. 

2  2  Fluid  Mechanics:  Fluid  flow  in  a  curved  thin 
film  is  also  controlled  by  the  pressure  along  the 
interface.  Using  the  definitions  for  q.4.  x0.  n0  and 
the  small  slope  assumption  given  above  leads  to 
d  3n0  n' 

d7(*K  +  tl)  =  if  (NcT»  ’  v>  (|5> 

Using  the  lubrication  approximation,  the  mass  flow 
rate  r  per  unit  width  in  a  slightly  tapered  liquid  film 
can  be  written  as 


S3  dP/ 
3v  dx 


(16) 


in  which  v  is  the  kinematic  viscosity. 

Using  Eqs.  (5.  15  and  16)  while  assuming  a  constant 
pressure  in  the  vapor  gives 


r  = 


83rio  q* 

- 2  (Ncq"  -  4) 

X0v  TT 


(17) 


from  which  follows  the  dimensionless  mass  flow  rate 
in  the  film 


The  transformed  differential  equation  is 

(^■-q3q’")  =  Nc  ( 1  -  3q "  -  q  3 )  (23) 

where  the  superscript  prime  now  refers  to 
differentiation  by  X. 

Applying  the  transformation.  Eq  (22).  to  the 
expression  for  the  total  mass  flow  rate  in  the  film.  Eq 
(18).  yields 


r  =  1'"  •  q>  ,:Jl 

where  the  superscript  prime  denotes  differentiation 
by  X.  If  the  meniscus  is  symmetric,  at  the  symmetry 

plane.  X  =  0.  the  mass  flow  rate  r  must  be  equal  to  zero 
with 


q  "  =  0  and  q'  =  0  ( 25  1 

Equation  (25)  provides  two  of  the  necessary  four 
initial  conditions.  The  values  of  q  and  q"  at  the  center 
plane  are  parameters. 

h  =  ho  X  =  0  (26) 
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Tl"  =  T)o 


X  =  0 


(27) 


These  initial  conditions  arc  different  from  the 
corresponding  conditions  used  in  previous  studies 
which  considered  a  meniscus  beginning  with  an 
adsorbed,  non-evaporating  thin  film  (i.e.,  r| (0)  =  I). 
First.  Gear's  method  was  applied  to  the  differential 
equation  (23)  and  boundary  conditions.  Eqs.  (25-27), 

for  a  given  value  of  Nc,  but  various  values  of  the 
initial  conditions  were  selected.  This  allows 
comparison  of  a  meniscus  with  a  small  adsorption 
effect  at  the  origin  (Case  A.  q(0)  =  4)  to  those 
previously  analyzed  with  a  large  adsorption  effect 
(Case  B.  rj (0)  =  1.02).  The  dimensionless  meniscus 
thickness  profile  and  dimensionless  evaporation  flux 
profile  arc  displayed  in  Figures  2  and  3.  These  figures 
suggest  that  a  cooling  device  using  a  meniscus 
between  two  feed  ports  (Figure  1)  with  a  minimum 
adsorption  effect  at  the  centerplanc  would  provide  a 
heat  flux  several  times  larger  than  a  device 
employing  a  meniscus  which  terminates  in  a  thin 
adsorbed  film.  In  the  case  at  hand  the  heat  flux  is 

roughly  8.8  times  larger.  This  result  reinforces  the 

view  that,  whereas  the  adsorption  effects  are 
necessary  to  keep  the  film  in  contact  with  the 
substrate  at  the  origin,  they  reduce  the  heal  flux. 
Further,  the  pressure  gradient  is  much  greater  in  the 
thinner  films.  In  both  cases  the  evaporation  rate  is 
fairly  constant.  However  in  case  A  »)'3  is  never  very 
large  whereas  in  case  B  starts  near  one  and 
decreases  to  near  zero.  Therefore  Equation  12  shows 
that  the  curvature  will  vary  weakly  in  case  A  but 
strongly  in  case  B.  In  case  B  it  increases  with 
increasing  X  and  then  decreases  slightly.  The 
increase  is  clearly  shown  by  Figure  2. 


Figure  7  Dimensionless  plot  ol  film  thickness 
versus  distance  lor  two  types  of 
meniscus  In  Case  A  p,,  =  4.(1  and  r|0  =0.15 
whereas  in  Case  B  Ho  -  1 .02  and  r|„  =  0  005 
In  each  case  \  -  I)  01 

A  comparison  was  also  developed  between 
solutions  having  different  values  of  Nc  but  the  same 
values  of  the  initial  conditions  (dimensionless).  These 
dimensional  solutions  are  shown  in  Figures  4-6  for 
the  meniscus  thickness,  curvature,  and  heat  flux 
profiles  for  the  example  system  of  n-octanc  on  silicon 
dioxide.  The  dimensional  pressure  at  the  ccnterplane 
is  the  same  in  each  case.  The  applicable  values  of  the 
physical  properties  and  groups  arc  given  in  Table  I. 


^nnn  4<oo  ejo  ^  12.00  16JJ0  ?o.oo 

Figure  3  Dimensionless  plot  of  evaporative  flux 
comparing  Case  A  with  Case  H 


0-0  X  ii^.m 

Figure  4  Dimensional  plot  of  him  thickness 

versus  distance  for  Case  I.  II.  and  III  i  see 
1  able  li.  In  cadi  case  no  =  4.1)  and  pu  - 


Figure  5 


Dimensional  plot  of  curvature  versus 
distance  for  Case  I.  II.  and  III 
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Figure  ft  Dimensional  plot  of  heal  flux  versus 
distance  for  Case  1.  II.  and  111 


Three  different  values  of  the  Hamaker  constant  were 
assumed  because  it  would  represent  a  major  variable 
in  real  systems.  The  value  of  the  Hamaker  constant 
has  been  calculated  and  measured  for  octane  on  SiC>2 
(Ingram,  1974;  Truong  and  Wayncr.  19871.  Due  to 
contamination  (e  g.,  an  adsorbed  water  film)  and/or  a 
slightly  microporous  surface  the  experimental  values 
have  been  found  to  be  the  lesser.  Therefore,  the 
theoretical  values  used  in  our  numerical  work  which 
lead  to  the  variation  in  Nc  covered  this  range.  The 
results  in  Figure  4  demonstrate  how  the  film  starts 
thinner  with  the  lower  value  of  the  Hamaker 
constant.  Figure  5  demonstrates  that  there  is  an 
increase  in  the  curvature  gradient  and,  therefore  the 
intcrfacial  mass  flux  with  a  decrease  in  the  value  of 
the  Hamaker  constant.  This  result  can  be  described 
using  Eq.  (8);  as  the  value  of  the  disjoining  pressure 
TABLE  1 


CASE 

i 

II 

III 

A.  1 

-10-22 

-10-23 

-10-24 

r.  n 

298 

298 

298 

AT,  °C 

0.00577 

0.00577 

0.00577 

Ns 

M.-T 
m  L 

5160  (10-7) 

5160  (10-7) 

5160  no-7) 

p.  kg/m1 

698.2 

698.2 

698.2 

IV  N/m2 

1.864  (103) 

1.864  (!03) 

1.864  (103) 

60.  IO-'°m 

27.296 

12.67 

5.881* 

Nc 

0.1 

0.464 

2.154 

N 

m 

21.3  (10-3) 

21.3  (10-3) 

21.3  (10-3) 

x0.m 

1.98  no-7) 

6.25  (10-8) 

1.98  (10-8) 

•Note:  5o  is  a  reference  thickness  only  and  does  not 
in  general  represent  the  exact  film  thickness. 


and/or  the  capillary  pressure  decreases,  the 
interfacial  mass  flux  increases.  The  pressure 
gradient  has  a  large  effect  and  an  increase  in  the 
value  of  Nc  signifies  an  increase  in  the  relative 
importance  of  capillarity.  Because  of  the  limitation 
on  the  range  of  validity  of  Eq.  (6)  lower  values  of  the 
temperature  difference  between  the  liquid  vapor 
interface  and  the  vapor  were  used.  These  lower 
values  increase  the  value  of  Nc.  reducing  the 
meniscus  thickness  in  the  transition  region  (the 
transition  region  bridges  the  adsorbed  thin  film 
dominated  by  disjoining  pressure  and  the  capillary 
meniscus  in  which  the  form  of  the  disjoining 
pressure  is  immaterial).  This  resulted  in  lower  values 
of  the  heat  flux. 

In  this  analysis  we  neglect  various  effects  and 
simplify  others.  The  temperature  of  the  vapor  liquid 
interface  is  taken  to  be  constant.  Thus  surface  stress 
due  to  a  variation  may  be  neglected.  The  accuracy  of 
the  assumption  may  be  checked  by  calculating  the 
temperature  drop  through  the  meniscus 


using  k  equal  0.14  W/mK.  Thus  a  film  thickness  may 
be  found  at  which  this  temperature  drop  is  10  percent 
of  the  temperature  difference  between  the  liquid  and 
the  vapor  phases.  The  results  for  cases  I.  II.  and  111 
arc  2760,  2410.  and  1380  Angstroms.  Therefore,  the 
film  resistance  is  apparently  negligible  for  a 
meaningful  portion  of  the  capillary  meniscus. 
However,  we  note  that  these  cases  were  not  chosen  to 
optimize  the  heat  flux,  but  to  add  insight  to  the 
phenomena.  Other  assumptions  include  lubrication 
flow,  and  the  form  of  disjoining  pressure.  Thus  the 
meniscus  must  be  fairly  flat  and  disjoining  effects 
must  be  negligible  where  it  is  thicker  than  200 
Angstroms.  Cases  II  and  III  satisfy  the  latter 
restriction  but  Case  I  is  borderline.  It  is  included  to 
illustrate  the  trend.  Another  model  is  needed  for 
thicker  films  (Truong  and  Wayncr.  1987). 

3.2  Additional  Analytical  Discussion 

The  model  of  Wayncr  ct  al.  (1976)  is  herein 
extended  to  include  capillary  effects.  The  results  may 
be  compared  in  the  adsorbed  thin  film  portion  of  the 
meniscus.  Here  both  models  have  analytical  solutions 
which  are  only  valid  if 

P/  •  P 

T|  -  1  and  — Ti — ~  -  -I  (29) 

1  'o 

The  isothermal  model  of  Rcnk  ct  al..  (1978)  also  has  an 
analytical  solution,  subject  to  the  same  conditions.  In 
each  of  these  three  cases  we  impose  the  boundary 
condition 

q  -»  1  as  (j  — >  -°o  (30) 

Our  result  for  the  present  model  was  developed 
with  Eq.  (21)  and  is 

T|  s  l  +  ac^/Nc,/2  +  bcW3  (31) 
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1 


3b(  1  -3NC )  c-5 


(32) 


Prrjaguin,  B  V  .  Nerpin.  S.V..  and  CTuraycv.  N.V..  1965, 
Ll'fcct  of  film  Transfer  upon  Evaporation  of  Liquids 
from  Capillaries.  Bull.  R.I.L.E.M..  vol.  29,  pp  93-98. 


r  =  -3V3  b(  1  -  Nc)  eW3  (33) 

where  a  and  b  are  arbitrary  constants 

The  model  of  Wayncr  el  al  (1976)  also  gives  (he 
same  results  as  Eqs.  (21),  (31),  (32).  and  (33)  with  the 
parameter  "a"  equal  to  zero.  Nc  set  equal  to  zero  in  Eqs. 
(32)  and  (33).  and  "b"  replaced  by  a  different 
arbitrary  constant.  B  The  non-evaporating  meniscus 
iRenk  el  al..  1978)  has  an  adsorbed  thin  film  as  well 
and  its  governing  equation  may  also  be  analyzed.  The 
governing  equation  is  given  by  setting  the  righthand 
side  of  Eq.  (21)  equal  to  zero  (this  follows  from  Eq.  (12) 
and  the  fact  that  there  is  no  evaporation  from  an 
isothermal  meniscus).  Our  rcsull  is  Eqs.  (31  >-( 33 )  with 
"b"  equal  zero  and  'a”  replaced  by  a  different 
arbitrary  constant  tr. 

In  the  earlier  model  of  Wayncr  ct  al.  ( 1976) 
capillary  effects  were  ignored.  Therefore  we  expect 
ilns  model  to  perform  best  if  the  parameter  Nc  is 
sin  j 1 1  In  (hat  case  (he  present  solution  for  q  w  ill 
always  be  equal  to  the  previous  solution  for  q  -  1 
sufficiently  small  because  the  lirst  exponential  term 
in  Eq.  (31)  may  be  neglected  in  favor  of  the  second 
term  if  \  is  sufficiently  negative  Of  course  the  more 
negative  \  is.  the  closer  q  is  to  I.  Under  these 
conditions,  the  solution  for  pressure  and  mass  flow 
are  almost  identical.  The  old  solution  is  merely 
corrected  by  a  term  proportional  to  N't.  According  to 
this  result,  the  effect  of  capillarity  is  found  to  slightly 
reduce  the  pressure  and  the  mass  flow. 

In  genera),  however,  the  previous  model  and 
present  model  agree  more  or  loss  Jcpcnding  on  the 
ratio  of  a/b  which  is  not  a  function  of  Nc.  The  details 
of  the  comparison  arc  beyond  the  scope  of  this  paper. 

4.  CONCLUSION 
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Herein,  we  have  demonstrated  the  effect  that 
pressure  and  temperature  can  have  on  the  heat 
transfer  characteristics  of  the  contact  line.  The 
el  feet  of  intcrfacial  forces  on  the  pressure  field  arc 

modeled  using  the  Humaker  constant  (6irA).  surface 
tension,  thickness,  and  curvature  At  a  given 
temperature,  the  pressure  field  affects  the  vapor 
pressure  and  therefore,  the  evaporation  rate.  In  the 
test  calculations,  the  evaporation  rate  was  increased 
by  decreasing  adsorption  Since  we  arc  evaluating 
the  effects  of  both  long  range  liquid  solid 
intcrmolccular  forces  and  liquid-vapor  intcrfacial 
tension  the  importance  of  these  results  naturally 
increases  as  the  system  becomes  smaller 
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